AFRL-HE-WP-TR-2005-0029 


United  States  Air  Force 
Air  Force  Materiel  Command 
Air  Force  Research  Laboratory 


Wave  Variable  Method  to  Control  Force-Reflecting 
Teleoperators  with  Time  Delays:  Generalizing  the  Wave 
Variable  Method  to  Multiple  Degree-of-Freedom 

Systems 

R. G.  Roberts 
C.A.  Moore 

Florida  A&M  University 
Tallahassee,  FL 

S.  Tosunoglu 

Florida  International  University 
Miami,  FL 

D.  W.  Repperger 
Air  Force  Research  Laboratory 
Human  Effectiveness  Directorate 
Wright-Patterson  AFB  OH 


March  2005 


Final  Report  for  October  2002  To  December  2004 


Approved  for  Public  Release; 
distribution  is  unlimited. 


20050526  053 


Human  Effectiveness  Directorate 
Warfighter  Interface  Division 
Collaborative  Interfaces  Branch 
2255  H  Street 

Wright-Patterson  AFB  OH  45433-7022 


STINFO  COPY 


NOTICES 


When  US  Government  drawings,  specifications,  or  other  data  are  used  for  any  purpose 
other  than  a  definitely  related  Government  procurement  operation,  the  Government 
thereby  incurs  no  responsibility  nor  any  obligation  whatsoever,  and  the  fact  that  the 
Government  may  have  formulated,  furnished,  or  in  any  way  supplied  the  said  drawings, 
specifications,  or  other  data,  is  not  to  be  regarded  by  implication  or  otherwise,  as  in  any 
maimer  licensing  the  holder  or  any  other  person  or  corporation,  or  conveying  any  rights 
or  permission  to  manufacture,  use,  or  sell  any  patented  invention  that  may  in  any  way  be 
related  thereto. 

Please  do  not  request  copies  of  this  report  from  the  Air  Force  Research  Laboratory. 
Additional  copies  may  be  purchased  from: 

National  Technical  Information  Service 
5285  Port  Royal  Road 
Springfield,  Virginia  22161 

Federal  Government  agencies  and  their  contractors  registered  with  the  Defense  Technical 
Information  Center  should  direct  requests  for  copies  of  this  report  to: 

Defense  Technical  Information  Center 
8725  John  J.  Kingman  Road,  Suite  0944 
Ft.  Belvoir,  Virginia  22060-6218 


TECHNICAL  REVIEW  AND  APPROVAL 
AFRL-HE-WP-TR-2005-0029 

This  report  has  been  reviewed  by  the  Office  of  Public  Affairs  (PA)  and  is  releasable  to 
the  National  Technical  Information  Service  (NTIS).  At  NTIS,  it  will  be  available  to  the 
general  public. 

This  technical  report  has  been  reviewed  and  is  approved  for  publication. 

FOR  THE  COMMANDER 

//Signed// 

MARIS  M.  VIKMANIS 

Chief,  Warfighter  Interface  Division 

Air  Force  Research  Laboratory 


REPORT  DOCUMENTATION  PAGE  OMBNa^olowa 


The  public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  Searching  existing  data  sources, 

gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of 
tformatlon,  including  suggestions  for  reducing  the  burden,  to  the  Department  of  Defense,  Executive  Services  and  Communications  Directorate  (0704-0188).  Respondents  should  be  aware 
that  notwithstanding  any  other  provision  of  law,  no  person  shed  be  subject  to  any  penalty  for  failing  to  comply  with  a  codaction  of  information  if  it  does  not  display  a  currently  valid  OMB 
control  number. 

PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ORGANIZATION. 


1.  REPORT  DATE  (DD-MM-YYYY)  2.  REPORT  TYPE 

01-03-2005  Final 


4.  TITLE  AND  SUBTITLE 

Wave  Variable  Method  to  Control  Force-Reflecting  Teleoperators  with  Time 
Delays:  Generalizing  the  Wave  Variable  Method  to  Multiple  Degree-of- 
Freedom  Systems 


6.  AUTHOR(S) 

R. G.  Roberts, A 

C. A.Moore.l 

S.  Tosunoglu, 

D. W.  Repperger ' 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

1  j 

Florida  A&M  University  Florida  International  University 

Wahnish  Way  &  Gamble  Street  CEAS,  IOSSS 

Tallahassee,  FL32307  West  Flagler  St,  EAS  2910 

Miami,  FL  33174 


9.  SPONSORING/MONITORING  AGENCY  NAME(S|  AND  ADDRESS(ES) 
Air  Force  Research  Laboratory, 3 
Human  Effectiveness  Directorate 
Warfighter  Interface  Division 
Air  Force  Materiel  Command 

Wright-Patterson  AFB,  OH  45433-7022 


12.  DISTRIBUTION/AVAILABIUTY  STATEMENT 

Approved  for  public  release;  distribution  is  unlimited. 


3.  DATES  COVERED  (From  -  To) 

Oct  02  to  Dec  04 


Sa.  CONTRACT  NUMBER 

F33615-02-1-6075 


Sb.  GRANT  NUMBER 

BAA-02-03-HEC-1 


6c.  PROGRAM  ELEMENT  NUMBER 

62202F 


64.  PROJECT  NUMBER 

7184 


6a.  TASK  NUMBER 

08 


Ef.  WORK  UNIT  NUMBER 

FL 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


10.  SPONSOR/MONITOR  S  ACRONYM(S) 

AERL/HECP 


11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 

AFRL-HE-WP-TR-2005-0029 


14.  ABSTRACT 

The  Wave  Variable  Method  is  used  to  stabilize  humans  interacting  with  teleoperation  tasks  apd  time  delays.  This  report  generalizes 
the  wave  variable  technique  tA  multiple  degree-of-ffeedom  systems.  A  new  class  of  scaling  matrices  are  obtained  which  has 
applicability  to  haptic  devices. 


16.  SUBJECT  TERMS 

Teleoperation,  Haptic  Interfaces,  Wave  Variable 


16.  SECURITY  CLASSIFICATION  OF:  | 

a.  REPORT 

U 

b.  ABSTRACT 

u 

c.  THIS  PAGE 

u 

17.  LIMITATION  OF 
ABSTRACT 


19a.  NAME  OF  RESPONSIBLE  PERSON 
Daniel  Repperger 


19b.  TELEPHONE  NUMBER  (Include  tree  code) 

937-255-8765 


Standard  Form  298  (Rev.  8/98) 

PrncrilMd  by  ANSI  Std.  Z39.1 8 


INSTRUCTIONS  FOR  COMPLETING  SF  298 


1.  REPORT  DATE.  Full  publication  date,  including 
day,  month,  if  available.  Must  cite  at  least  the  year 
and  be  Year  2000  compliant,  e.g.  30-06-1998; 
xx-06-1998;  xx-xx-1998. 

2.  REPORT  TYPE.  State  the  type  of  report,  such  as 
final,  technical,  interim,  memorandum,  master’s 
thesis,  progress,  quarterly,  research,  special,  group 
study,  etc. 

3.  DATES  COVERED.  Indicate  the  time  during 
which  the  work  was  performed  and  the  report  was 
written,  e.g.,  Jun  1 997  -  Jun  1 998;  1-10  Jun  1 996; 
May -Nov  1998;  Nov  1998. 

4.  TITLE.  Enter  title  and  subtitle  with  volume 
number  and  part  number,  if  applicable.  On  classified 
documents,  enter  the  title  classification  in 
parentheses. 

5a.  CONTRACT  NUMBER.  Enter  all  contract 
numbers  as  they  appear  in  the  report,  e.g. 
F33615-86-C-5169. 

5b.  GRANT  NUMBER.  Enter  all  grant  numbers  as 
they  appear  in  the  report,  e.g.  AFOSR-82-1234. 

5c.  PROGRAM  ELEMENT  NUMBER.  Enter  all 
program  element  numbers  as  they  appear  in  the 
report,  e.g.  61 101  A. 

5d.  PROJECT  NUMBER.  Enter  all  project  numbers 
as  they  appear  in  the  report,  e.g.  1F665702D1257; 
ILIR. 

9 

5e.  TASK  NUMBER.  Enter  all  task  rfumbers  as  they 
appear  in  the  report,  e.g.  05;  RF0330201;  T41 12. 

'  * 

5f.  WORK  UNIT  NUMBER.  Enter  ail  work  unit 
numbers  as  they  appear  in  the  report,  e.g.  001; 
AFAPL304801 05. 

6.  AUTHOR(S).  Enter  name(s)  of  person(s) 
responsible  for  writing  the  report,  performing  the 
research,  or  credited  with  the  content  of  the  report. 
The  form  of  entry  is  the  last  name,  first  name,  middle 
initial,  and  additional  qualifiers  separated  by  commas, 
e.g.  Smith,  Richard,  J,  Jr. 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND 
ADDRESS(ES).  Self-explanatory. 


8.  PERFORMING  ORGANIZATION  REPORT  NUMBER. 
Enter  all  unique  alphanumeric  report  numbers  assigned 
by  the  performing  organization,  e.g.  BRL-1234; 
AFWL-TR-85-4017-Vol-21-PT-2. 

9.  SPONSORING/MONITORING  AGENCY  NAME(S) 
AND  ADDRESS(ES).  Enter  the  name  and  address  of  the 
organization(s)  financially  responsible  for  and 
monitoring  the  work. 

10.  SPONSOR/MONITOR'S  ACRONYM(S).  Enter,  if 
available,  e.g.  BRL,  ARDEC,  NADC. 

11.  SPONSOR/MONITOR’S  REPORT  NUMBER(S). 

Enter  report  number  as  assigned  by  the  sponsoring/ 
monitoring  agency,  if. available,  e.g.  BRL-TR-829;  -215. 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT.  Use 
agency-mandated  availability  statements  to  indicate  the 
public  availability  or  distribution  limitations  of  the 
report.  If  additional  limitations/  restrictions  or  special 
markings  are  indicated,  follow  agency  authorization 
procedures,  e.g.  RD/FRD,  PROPIN,  ITAR,  etc.  Include 
copyright  information. 

13.  SUPPLEMENTARY  NOTES.  Enter  information  not 
included  elsewhere  such  as:  prepared  in  cooperation 
with;  translation  of;  report  supersedes;  old  edition 
number,  etc. 

14.  ABSTRACT.  A  brief  (approximately  200  words) 
factual  summary  of  the  most  significant  information. 

15.  SUBJECT  TERMS.  Key  words  or  phrases 
identifying  major  concepts  in  the  report. 

16.  SECURITY  CLASSIFICATION.  Enter  security 
classification  in  accordance  with  security  classification 
regulations,  e.g.  U,  C,  S,  etc.  If  this  form  contains 
classified  information,  stamp  classification  level  on  the 
top  and  bottom  of  this  page. 

17.  LIMITATION  OF  ABSTRACT.  This  block  must  be 
completed  to  assign  a  distribution  limitation  to  the 
abstract.  Enter  UU  (Unclassified  Unlimited)  or  SAR 
(Same  as  Report).  An  entry  in  this  block  is  necessary  if 
the  abstract  is  to  be  limited. 


Standard  Form  298  Back  (Rev.  8/98) 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


11 


Acknowledgements 


This  work  was  partially  funded  by  the  U.S.  Air  Force  Research  Laboratory 
(AFRL/HECP)  Human  Effectiveness  Directorate,  Integration  and  Operations  Divisions 
(AFRL/HEO);  Award  NO:  F33615-02-1-6075.  The  funding  was  provided  to  Florida 
International  University  (FIU)  in  Miami,  and  Florida  A&M  University  (FAMU)  in 
Tallahassee,  where  both  teams  collaborated  during  the  investigation  in  conjunction  with 
the  technical  and  managerial  supervision  of  Dr.  Daniel  W.  Repperger.  FIU  and  FAMU 
teams  acknowledge  this  funding  without  which  this  work  would  not  be  possible. 


iii 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


IV 


GENERALIZING  THE  WAVE  VARIABLE  METHOD  TO  MULTIPLE 
DEGREE-OF-FREEDOM  SYSTEMS 

I.  INTRODUCTION 

One  of  the  tasks  of  the  project  was  to  do  a  literature  survey.  Various  articles  were  examined, 
including  references  [l]-[6].  We  examined  in  detail  some  of  the  work  done  by  Wayne  Book  and 
Sahgir  Munir  at  Georgia  Institute  of  Technology.  In  their  work,  they  generalized  the  wave 
variable  method  to  multiple  degree-of-freedom  teleoperation  systems  by  replacing  the  damping 
coefficient  in  the  standard  wave  variable  method  with  a  set  of  scaling  matrices.  In  the  following, 
we  will  derive  a  more  general  set  of  scaling  matrices.  In  the  next  section,  we  review  the  single 
degree-of-freedom  and  multiple  degree-of-freedom  systems.  In  Section  III,  we  derive  a  larger 
family  of  scaling  matrices  and  verify  that  an  important  subset  of  scaling  matrices  results  in  stable 
operation  of  a  multiple  degree-of-freedom  teleoperation  system  with  a  fixed  but  unknown  time 
delay.  We  then  determine  the  relationship  between  the  extended  family  of  scaling  matrices  and 
the  family  proposed  by  Munir  and  Book.  Section  IV  contains  a  simulation  to  illustrate  the 
concepts  and  Section  V  contains  the  conclusions. 

n.  THE  WAVE  VARIABLE  METHOD  FOR  MULTIPLE  DEGREE-OF-FREEDOM 

SYSTEMS 


The  wave  transformation  relations  are  given  by 

uAt)  =  um(t-T) 
v«(0  =  v,(r-r). 

The  wave  transformations  for  the  left  wave  junction  are  given  by 

4>0Jt)  +  Tm(t) 


Vm(0  = 


42b 

bem(t) 


(i) 


(2) 


*■«(?) 


42 b 


and  that  for  the  right  wave  junction  are  given  by 

Ww(|)+UO 


«,(0=- 


v,(0=- 


42b 


(3) 


42b 


Although  the  strictly  positive  parameter  b  can  be  chosen  arbitrarily,  it  defines  a  characteristic 
impedance  associated  with  the  wave  variables  and  directly  affects  the  system  behavior  [5]. 

Equations  (2)  and  (3)  are  for  1  DOF  systems.  To  implement  the  wave  variable  method  on  a 
system  that  has  more  than  one  degree  of  freedom,  the  equations  for  the  transforms  must  be 
generalized.  Niemeyer  and  Slotine  [5]  suggest  making  b  a  positive  definite  matrix.  Munir  and 
Book  [2]-[4]  recommend  a  more  general  formulation  by  writing  the  transformation  equations 
from  before  as 

vs(t)  =  Cjs(t)-DwTlAt) 
and 


(4) 


I 


(5) 


v„(0  = 


«tW=AA(')  +  V#„,(0 

where  /lu.,  Bn,  CK,  and  D„.  are  nXn  scaling  matrices  and  n  is  the  number  of  degrees  of 

freedom  of  the  teleoperation  system.  These  matrices  cannot  be  chosen  arbitrarily;  it  is  necessary 
to  determine  conditions  for  the  scaling  matrices  to  guarantee  passivity  [6].  To  accomplish  this  we 
will  define  the  power-flow  at  each  side  to  be 

K,  Tn,  =  l/2  «>,„  -  Y2  vl  vm  (6) 

for  the  master  side  and 


dy,,,!  =  + Yi  vi y. 


(7) 


for  the  slave  side.  Substituting  equations  (4)  and  (5)  into  equation  (6)  or  (7),  expanding,  and 
matching  matrix  coefficients  yields  the  requirements 


AlAw=ClCw 
BTWBW  =  DtwDw 


(8) 


and  also  that 


I=/2(2AlBw+2ClDw). 


(9) 


Munir  and  Book  [2]-[4]  derive  conditions  on  Aw  and  Biv  to  ensure  that  (8)  and  (9)  are  satisfied. 
In  particular,  they  note  that  the  scaling  matrices  must  be  nonsingular  and  consider  the  special  case 


Aw  —  Cw 
bk  =  dw 


(10) 


so  that  equation  (8)  is  satisfied.  Using  this  relationship,  equation  (9)  reduces  to 

I  =  2AlBw.  (11) 

Munir  and  Book  then  restrict  Aw  to  be  symmetric  but  not  necessarily  positive  definite  and  prove 

that  the  resulting  family  of  scaling  matrices  results  in  a  stable  system  by  showing  that  passivity  is 
ensured  using  the  norm  of  the  scattering  matrix. 

Although  this  is  more  general  than  what  was  proposed  in  [5],  because  of  the  specific  choices 
made,  only  a  restricted  class  of  scaling  matrix  is  determined.  In  this  article,  we  will  extend  the 
family  of  scaling  matrices  that  result  in  stability  and  discuss  the  significance  of  this  extension  in 
terms  of  the  wave  variables  themselves. 


III.  DERIVATION  OF  A  LARGER  FAMILY  OF  SCALING  MATRICES 


To  extend  the  family  of  scaling  matrices  originally  proposed  by  Munir  and  Book,  we  will  first 
derive  the  whole  family  of  matrices  satisfying  (8)  and  (9).  From  equation  (8)  it  is  clear  that  Al 

and  Cl  have  the  same  column  space  and  that  BK  and  £)„.  have  the  same  row  space.  This 

observation  along  with  equation  (9)  implies  that  all  four  scaling  matrices  must  be  nonsingular. 
Since  all  square  root  decompositions  of  nonsingular  square  matrices  such  as  those  given  in  (8)  are 
related  by  pre-multiplication  by  an  orthogonal  matrix,  it  follows  that 

=  (7|  Aw 

=  Q2Bw 

where  Qt  and  Q2  are  nXn  orthogonal  matrices.  Now  this  is  not  enough  to  guarantee  that  (9)  is 
satisfied.  We  will  say  that  Alt,  and  Bw  are  compatible  with  respect  to  (8)  and  (9)  if  there  are 
orthogonal  matrices  g,  and  Q2  so  that  A,„,  Bw,  CK  =Q,Air,  and  Dn  =Q2Bu  satisfy  (8)  and  (9). 


(12) 
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Substituting  (12)  into  (9)  and  performing  the  required  manipulations  gives  the  following 
necessary  and  sufficient  condition  for  Aw  and  BK  to  be  compatible  in  this  sense: 

A  BT  +B  Ar  =1  .  (13) 

w  w  w  w  v  ' 

With  some  algebraic  manipulation,  it  can  then  be  shown  that  the  characterizing  condition  for  Aw , 


Bw,  CK ,  and  Dw  to  satisfy  (8)  and  (9)  is  that  the  following  hold: 


1.  Aw  is  nonsingular. 

2.  Bw  =  j/^(I  +  S)A~T  where  S  is  any  nxn  skew-symmetric  matrix. 

3.  CB.  =  £M„.  where  Q  is  any  nxn  orthogonal  matrix. 

4-  Dw=y2Q(I-S)Af. 

Note  that  these  four  conditions  guarantee  that  all  four  matrices  are  nonsingular. 

Next,  we  identify  a  subfamily  of  scaling  matrices  that  result  in  stable  operation.  To  do  this,  we 
need  to  determine  the  scattering  matrix.  The  transfer  function  for  the  input-output  relationship 
across  the  communication  link  can  be  shown  to  be  [2] 


GJs) 


1  +  e~2xT  " 
2e~sT  , 


-  2e~sT  f  .1 


(1 


1  +  e 


-2  sT 


1  +  e 


1  +  e~‘ 

■a:b. 


-2xT 


(14) 


GK(jco)  = 


(15) 


Setting  s  =  j(0  and  simplifying  gives 

B~}AW  ta.nh(  jcoT)  -sech (  jcoT)I  1 

sech(;<yr)7  A~'BIV  tanh(yW)J 
Passivity  can  then  be  demonstrated  using  the  scattering  matrix 

S{j(o)  =  [G„  (jco)  -  I][GuXj(0)  +  /]-’.  (16) 

While  calculations  involving  the  inversion  of  partitioned  matrices  are  generally  difficult  to  do  in 
closed  form,  (15)  has  enough  mathematical  structure  to  allow  this.  It  can  be  shown  that  the 
expression  for  S(jco)  is 

~[M 2  -I]sinh(  jcoT)  -2  M  1 

2  M  -[M1  -/]sinh(yW)J 

where  M  -  B~?AW  and  where  D(joJT)  is  given  in  block  diagonal  form  as 

D(jcoT)  =  diag(D0  (jcoT),  D0  (jcoT))  ( 1 8) 

with 


D(jcoT) 


(17) 


Dtl(jwT)  =  [2M  cosh(jcoT)  +  (M2  +  I)sinh(jcoT)Y'  ■  (19) 

Since  the  calculation  of  D0(jcoT)  requires  a  matrix  inversion  of 

Kx  cos(ftiT)  +  jK2  sin(cJT)  ^20) 

=  2Mcos(coT)  +  j[M 1  +  /]sin(ft>T), 

it  must  be  shown  that  (20)  is  nonsingular  for  all  co.  This  verification  for  the  whole  class  of 
scaling  matrices  satisfying  (8)  and  (9)  requires  more  work  than  the  subfamily  proposed  by  Munir 
and  Book  [2]-[4],  The  approach  we  have  taken  is  to  show  that  (13)  imposes  restrictions  on  the 
eigenvalues  of  B~lAw  that  preclude  (20)  from  being  singular.  This  was  done  by  showing  that  if 

(20)  is  singular  then  B~'AW  necessarily  has  an  eigenvalue  with  zero  real  part.  However,  it  can  be 
shown  that  for  our  family  of  scaling  matrices  the  eigenvalues  of  B~'AU,  necessarily  have  positive 
real  parts. 
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Now  if  (20)  is  singular  then  so  is  its  product  with  its  complex  conjugate: 

[AT,  cos (coT)  +  jKj  sin(&)7)][  A",  cos(®7')  -  jK2  sin(ftjT)]  . 

=  K?  cos2  (afT)  +  Kl  sin2  (a)T)  +  j[K2K,  -  KlK2]sin(coT)cos(coT) 

=  sin2  (0)T)(B;J  A„Y  +  [2 +  2  cos1  (a)T)](B;' A  J2  +sin \o)T)l. 

Note  that  the  resulting  imaginary  part  is  zero  precisely  because  K]  =  2 M  and  K2  =  M 1  + 1 
commute;  otherwise  (21)  would  be  complex.  We  will  show  that  this  matrix  is  nonsingular  and  in 
the  process  prove  that  (20)  is  also  nonsingular.  If  sin(<yT)  =  0  then  (21)  becomes  4 (B~lAw)2 
which  is  nonsingular.  Suppose  that  sin(<s/T)^0.  In  order  for  (21)  to  be  singular  at  least  one  of 
the  four  roots  of 

p(A)  =  sm2  (coT)  A4  +  [2  + 2cos2  (coT)]A2 +s\n2  (coT)  (22) 

must  be  an  eigenvalue  of  B~'  Aw .  The  roots  of  this  polynomial  (which  is  a  quadratic  in  /l2)  are 
given  by 

A  = 

Since  (2 +  2  cos2  (coT))2  dominates  4sin4(ruT),  the  inner  square  root  is  real.  Furthermore, 
-2-2cos 2(q)T)  dominates  the  inner  square  root  so  the  number  inside  the  outer  square  root  is 
real  and  nonpositive.  Therefore,  the  roots  of  p(A)  are  purely  complex  (i.e.,  have  zero  real  part). 

Now  the  Aw  and  Bw  matrices  satisfying  (8)  and  (9)  also  satisfy  (13).  Because  Aw  and  Bw  are 
invertible,  (13)  can  be  written  as 

b-Jaw  +  aIb;J  =b:!b;J.  (24) 

This  says  that  the  symmetric  part  of  B~'AW  is  positive  definite  implying  that  the  real  parts  of  the 
eigenvalues  of  B~'AW  are  strictly  positive  so  that  none  of  the  eigenvalues  of  B~'  Aw  are  roots  of 
the  polynomial  p(A) .  Therefore  (20)  is  nonsingular  implying  that  Du(ja>)  is  weH  defined. 

Munir  [2]  demonstrates  that  the  system  is  stable  for  their  family  of  scaling  matrices  by  showing 
that  the  norm  of  the  scattering  matrix 

||5||  =  supAI/2(5*(7ty)5(jiy))  (25) 

eo 

is  equal  to  1  where  S'(ja) )  denotes  the  complex  conjugate  transpose  of  S(  jco) .  This  was  done 
by  showing  that 

S\jo))S(jco)  =  /  .  (26) 

However,  for  this  calculation  to  work  out,  it  must  be  assumed  that  M  =  B~'AW  is  symmetric, 
which  is  true  for  the  scaling  matrices  in  [2]-[4]  by  the  symmetry  of  Au,  along  with  equation  (1 1). 

If  we  constrain  the  family  of  matrices  derived  earlier  that  satisfy  equations  (8)  and  (9)  to  have 
the  property  that  B~'  Aw  is  symmetric,  the  same  calculations  used  by  Munir  [2]  to  ensure  passivity 

for  his  family  of  scaling  matrices  would  apply  to  give  (26)  for  our  constrained  family  of  matrices. 
This  will  of  course  reduce  the  size  of  our  family  of  scaling  matrices.  To  see  this,  note  that 

b;'ah  =  [1/2(I  +  S)A;T  r‘  au.  =  2  Al  (/ + sy'  a,„  ■  (27) 

Now  B~'AW  is  symmetric  if  and  only  if  its  inverse  is.  This  inverse  is  given  by 

=y2K'd +S)K?  =y2K'AT  +y2K'sA;!  ■  (28) 
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The  matrix  V  A~'A~J  is  symmetric  while  the  matrix  \/  A~'SA~J  is  skew-symmetric  so  in  order 

for  (28)  to  be  symmetric  we  have  that  S  =  0 .  This  is  the  necessary  and  sufficient  condition  for 
B~'AW  to  be  symmetric  for  the  class  of  matrices  satisfying  (8)  and  (9).  Thus  the  following  family 
of  scaling  matrices  can  be  used  for  applying  the  wave  variable  method  to  a  multiple  degree-of- 
freedom  system: 

1 .  Aw  is  nonsingular. 

2.  Bn=/2A^. 

3.  C„.  =  QAK  where  Q  is  any  nxn  orthogonal  matrix. 

4.  DK,=y2QA;T. 

Note  that  setting  Aw  =  and  Q  =  1  for  the  scalar  case  results  in  exactly  the  same  solution  as 


(2)  and  (3). 


Since  choosing  a  set  of  scaling  matrices  requires  the  selection  of  an  nxn  nonsingular  matrix  Aw 
and  an  nxn  orthogonal  matrix  Q,  there  are  a  total  of  n2  +n(n-\)/2  =  3n2 /2-n/2  degrees  of 
freedom  in  choosing  the  scaling  matrices.  However,  because  the  collection  of  scaling  matrices 
presented  in  [2]-[4]  are  uniquely  defined  by  the  selection  of  a  single  symmetric  matrix  A, the 

family  of  scaling  matrices  in  [2]-[4]  is  only  n(n  + 1)/2  -dimensional. 

It  is  then  natural  to  ask  how  this  extension  of  the  scaling  matrices  affects  the  wave  variables. 
To  see  this,  we  first  consider  the  effect  of  Q .  The  v  wave  variables  are  given  in  equations  (4)  and 
(5).  Substituting  in  the  expressions  for  CH.  and  Dn.  yields 

vJt)  =  QAjm(t)-y2QA;TTjt)  (29) 

=Q[Aj,„(t)-y2A:TTjt)i 

which  clearly  demonstrates  that  Q  merely  applies  an  orthogonal  transformation  to  the  v -variable, 
i.e.,  it  will  merely  rotate  and/or  reflect  the  v -variable.  The  same  holds  for  vt(f) .  This  will  clearly 

have  no  effect  on  the  power  flow  equations  (6)  and  (7). 

Next,  consider  the  effect  of  allowing  the  matrix  Aw  to  be  nonsymmetric.  There  is  a  well- 

known  result  in  matrix  theory  called  the  polar  decomposition  that  states  that  any  square  matrix 
M  can  be  written  as  a  product  M  =  PUT  where  P  is  a  symmetric  positive  definite  matrix  and 
U  is  an  orthogonal  matrix.  Setting  M  =  Al  gives  that  Aw  -  UPW  for  some  suitable  orthogonal 


matrix  U  and  symmetric  positive  definite  matrix  Pw .  The  u  wave  variable  then  becomes 

«„  (0  =  QPjm  (0  +  X {QP" yT  T" {t)  (30) 

=Q[pA(t)+y2pfTjt)] 

so  that  the  m  -variable  is  merely  rotated  and/or  reflected.  A  similar  statement  holds  for  the  v- 
variable. 


IV.  SIMULATION  RESULTS 


The  performance  of  the  wave  variable  method  was  tested  using  a  2  degree-of-freedom  linear 
model.  The  equations  of  motion  for  both  the  master  and  slave  manipulators  were  given  by 

Tri  =  j,e+B,e  (3i) 
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where  rpil  is  the  input  torque,  J%  is  the  desired  2x2  constant  symmetric  positive  definite  inertia 
matrix,  and  Bs  is  the  damping  matrix  with  the  same  qualities  as  7  .  With  the  wave  variable 


parameters  and  the  equations  of  motion  necessary  to  complete  the  system,  we  can  simulate  a  2 
degree-of-freedom  bilateral  teleoperation  system. 

As  can  be  seen  in  Fig.  2,  the  behavior  of  the  linear  2  degree-of-freedom  slave  is  similar  to  the 
behavior  of  the  master.  In  this  simulation,  the  input  torques  were  3  and  -4  Nm,  respectively.  The 
slave  manipulator  was  able  to  match  the  performance  of  the  master  manipulator  with  respect  to 
torque,  position,  and  velocity.  Also,  in  less  than  4  seconds  the  output  torques  had  reached  their 
steady  state  values.  Due  to  the  nature  of  the  wave  variable  method,  it  can  be  shown  that  the 
system  will  stabilize  for  any  amount  of  constant  time  delay. 


Fig.  2:  The  response  of  a  2-DOF  system  with  a  2T  =  1  sec  total  time  delay. 


Fig.  3:  A  plot  of  the  wave  variable  v,  for  the  2-DOF  teleoperation  system  for  two  different  q 
matrices. 

Fig.  2  represents  the  simulation  results  for  a  given  fixed  Aw  where  Q  is  allowed  to  be  any 
orthogonal  matrix.  No  matter  which  Q  is  used,  the  overall  system  output  will  not  change.  What 
do  change  are  the  wave  variables  that  are  sent  across  the  communication  line.  Fig.  3  shows  a  plot 
of  the  wave  variable  v,  for  the  cases  when  Q  is  the  identity  matrix  and  when  it  is  a  rotation 
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matrix  corresponding  to  a  rotation  of  45  degrees.  The  figure  clearly  shows  that  v,  will  be  rotated 
by  the  same  amount  that  CK  and  Dw  are.  Also,  when  looking  at  Pu  =UTA„,  which  is  a  rotated 
version  of  Aw,  the  same  statement  from  before  can  be  made  again.  This  time  however  both  um 
and  vt  will  be  rotated,  and  again  there  will  be  no  change  to  the  system  output. 

V.  CONCLUSIONS 

In  summary,  we  have  presented  the  derivation  of  an  extension  of  the  wave  variable  method  to 
multiple  degree-of-freedom  systems.  We  have  shown  that  the  collection  of  scaling  matrices 
determined  in  [2]-[4]  to  preserve  passivity  is  not  complete  and  have  determined  a  larger  family  of 
feasible  scaling  matrices.  We  have  also  shown  how  the  new  scaling  matrices  relate  to  those 
proposed  in  [2]-[4]. 
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